Broadening of the Beresinkii-Kosterlitz-Thouless superconducting transition by 

inhomogeneity and finite-size effects 
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We discuss the crucial role played by finite-size effects and inhomogeneity on the Beresinkii- 
Kosterlitz-Thouless (BKT) transition in two-dimensional superconductors. In particular, we focus 
on the temperature dependence of the resistivity, that is dominated by superconducting fluctuations 
above the BKT transition temperature Tbkt and by inhomogeneity below it. By means of a 
renormalization-group approach we establish a direct correspondence between the parameter values 
used to describe the BKT fluctuation regime and the distance between Tbkt and the mean-field 
Ginzburg-Landau transition temperature. Below Tbkt a resistive tail arises due to finite-size effects 
and inhomogeneity, that reflects also on the temperature dependence of the superfluid density. We 
apply our results to recent experimental data in superconducting LaAlOs/SrTiOs heterostructures, 
and we extract several informations on the microscopic properties of the system from our BKT fitting 
parameters. Finally, we compare our approach to recent data analysis presented in the literature, 
where the physical meaning of the parameter values in the BKT formulas has been often overlooked. 

PACS numbers: 74.20.-z, 74.25.Fy, 74.78.Fk 



I. INTRODUCTION 

In the last years a renewed interest emerged in the 
superconducting transition in two-dimensional (2D) sys- 
tems, prompted by the experimental achievement of high- 
quality ultra-thin films of superconducting materials. To 
this category belong both few-unit-cell thick films of 
layered cuprate superconductors 1-3 and the nanometer- 
thick layers of superconducting electron systems formed 
at the interface between insulating oxides in artificial 
LaA10 3 /SrTi0 3 heterostructures 4 ' 5 . At the same time, 
the experimental progresses made in the last decade of 
intense investigation in high-temperature superconduc- 
tors prompted additional measurements in thin films of 
conventional superconductors by means of different tech- 
niques or higher resolution than the ones available in 
the past. Typical examples are provided by the finite- 
frequency study of the optical magneto-conductivity in 
films of InOj; 6 , Nernst-effect measurements in amorphous 
films of NbSi 7 and scanning tunneling microscopy in TiN 
films 8,9 . 

Due to the 2D nature of these systems, the super- 
conducting (SC) transition is expected to belong to 
the Beresinskii-Kosterlitz-Thouless (BKT) universality 
class 10-12 , where the gauge symmetry is unbroken in the 
SC state but the system has a finite superfluid density, 
which is destroyed at Tbkt by proliferation of vortex- 
antivortex phase fluctuations. As it has been discussed 
at length in the past literature, the BKT transition has 
in principle very specific signatures 13,14 . For example, 
by approaching the transition from below, the superfluid 
density n s is expected to go to zero discontinuously at 
the BKT temperature Tbkt, with an "universal" relation 
between ii s (Tbkt) and Tbkt itself 13-15 . Approaching 



instead the transition from above one has in principle 
the possibility to identify the BKT transition from the 
temperature dependence of the superconducting fluctua- 
tions. Indeed, in 2D the temperature dependence of sev- 
eral physical quantities (like the paraconductivity or the 
diamagnetism) is encoded in the temperature dependence 
of the superconducting correlation length £(T), that di- 
verges exponentially at Tbkt, in contrast to the power- 
law expected within Ginzburg-Landau (GL) theory 16 . A 
possible interpolation scheme between standard GL fluc- 
tuations and BKT phase fluctuations of the SC order 
parameter was proposed long ago in a seminal paper by 
Halperin and Nelson 17 (HN). 

In the attempts made in the past to find out exper- 
imental signatures of the BKT transition in thin films 
of conventional superconductors 18-20 it turned out that 
BKT fluctuations are usually restricted to a small tem- 
perature regime near the GL transition temperature T c . 
If the energy range T c — Tbkt is extremely small, most of 
the fluctuation regime is dominated by GL fluctuations, 
and the exponential signatures of BKT fluctuations can 
be hardly detected. The predominance of the GL fluc- 
tuation regime has been confirmed also by more recent 
measurements of Nernst effect in NbSi films 7 and zero- 
bias tunneling conductance in TiN films 9 . Analogously, 
the expected universal jump of n s (TBKT) at Tbkt due 
to vortex proliferation can be overscreened by the simul- 
taneous fast decrease of n s (T) due to quasiparticle exci- 
tations near T c 20 . 

An additional effect that can mask the occur- 
rence of BKT transition is the intrinsic inhomogene- 
ity of the sample. For example, as it has been dis- 
cussed recently in the context of thin films of high- 
temperature superconductors 21 , the spatial inhomogene- 
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ity can broaden considerably the universal jump of the 
superfluid density, leading to a smooth downturn of the 
superfluid density instead of the sharp one expected in 
ultra-thin samples. Such an intrinsic mesoscopic inho- 
mogeneity has been revealed by scanning tunneling spec- 
troscopy in cuprate superconductors 22 , and recently also 
in films of conventional supercondutcors 8 . This indicates 
that inhomogeneity is a crucial ingredient to several su- 
perconducting systems in the presence of disorder, as sug- 
gested also by recent numerical simulations 23 . Finally, 
one must account also for finite-size effects (at the scale 
of the sample dimensions or even smaller, in the inhomo- 
geneous case) , that are expected to cut off the divergence 
of the correlation length at Tbkt- 

In this work we aim to address the role of inhomogene- 
ity and finite-size effects in the BKT transition, provid- 
ing a general scheme for analyzing paraconductivity mea- 
surements in quasi-2D superconductors. We first match 
the HN interpolation scheme with a detailed analysis of 
the renormalization group (RG) equations for the BKT 
transition far from the critical region where an analyti- 
cal solution is available 12 , taking into account finite-size 
effects. This analysis allows us to relate the parame- 
ters of the standard BKT correlation-length expression 
to microscopic quantities, reducing considerably the de- 
grees of freedom in the fitting procedure of the resistivity 
data. In particular, we show that once fixed the differ- 
ence T c — Tbkt and the value of the vortex-core energy, 
whose role has been recently discussed in the context of 
the physics of high-temperature superconductors 21,24-27 , 
the behavior of the resistivity from Tkt up to tem- 
peratures far above T c is uniquely determined. In this 
way we establish a consistency check that the fitting pa- 
rameters must satisfy when the standard BKT approx- 
imated formulas are used, a fact that has been often 
overlooked in the literature. Building on such an anal- 
ysis we can also take into account the role of inhomo- 
geneities, which as we will show are crucial to under- 
stand the experimental situation. Indeed even relatively 
small inhomogeneities can provide a significant enhance- 
ment of the resistive response below the Tbkt transi- 
tion, that is simultaneously reflected in the temperature 
dependence of the superfluid density across the transi- 
tion itself. A paradigmatic example of application of our 
analysis is provided by recent measurements in supercon- 
ducting heterostructures 4,5 . As we shall see, a correct 
treatment of finite-size effects and inhomogeneity allows 
us to reproduce with great accuracy the available exper- 
imental results, and to estimate the superfluid density in 
these unconventional systems. 

We note that although the BKT transition has been al- 
ready invoked to discuss the physics of such heterostruc- 
tures, a very different approach was taken so far in the 
literature 4,5,28 , relying essentially on the idea that the GL 
temperature T c is far larger than Tbkt (that can have 
eventually a different qualitative character 4 ), so that the 
whole fluctuation regime can be described within BKT 
theory. A schematic view of the difference between our 



analysis and the previous approach is shown in Fig. 1. 
The point of view taken in Refs. [4,5,28] leads of course 
to a very different identification of the various tempera- 
ture regimes, and very different physical parameters for 
the underlying BKT theory. In our view these previous 
analyses suffer from several problems, that we will dis- 
cuss explicitly in this paper. In particular they do not 
correctly identify the physical parameters of the system 
and thus overlook several interesting consequences that 
one can extract from the accurate experimental measure- 
ments performed in these new materials. 
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FIG. 1: (Color online) Comparison between the approach 
discussed in our paper and the one presented so far in the 
literature 4,5,28 as far as the resistivity data in superconduct- 
ing heterostructures are concerned. The experimental data 
for the resistivity R normalized to the normal-state value Rn 
are taken from Ref. 5, and are the same as the ones showed 
in Fig. 7 and Fig. 9 below. On the left panel we summarize 
our approach: most of the fluctuation regime above Tbkt is 
dominated by GL fluctuations, while BKT fluctuations are re- 
stricted to a narrow range of temperatures near the transition 
temperature Tbkt that would be observable in the homoge- 
neous system. However, inhomogeneity leads to a consider- 
able tail in the resistivity, that remains finite below Tbkt- 
On the right panel we summarize the point of view used in 
Ref. 4,5,28: the whole range of temperatures above Tbkt 
is dominated by SC fluctuations having BKT character, and 
finite-size effects are responsible for the resistive tail below 
Tbkt- 

The structure of the paper is the following. In Sec. II 
we review the standard description of fluctuation conduc- 
tivity in 2D superconductors, that establishes the corre- 
spondence between resistivity and the fluctuation corre- 
lation length. In Sec. Ill we analyze systematically the 
behavior of the correlation length within BKT theory, by 
means of a RG approach. This allows us to fully identify 
the parameters in the HN formula interpolating between 
the BKT transition and the standard Gaussian fluctua- 
tions. Sec. IV is devoted to the discussion of the role of 
inhomogeneity, trough the direct application of our inter- 
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polating BKT-to-GL scheme to the resistive transition. 
We apply our scheme to identify the relevant parame- 
ters in superconducting heterostructures. In Sec. V we 
clarify the differences between our approach and previous 
theoretical reports. The reader interested only to the is- 
sue of analyzing the experimental data can refer directly 
to the Sections IV and V. Finally, Sec. VI contains the 
concluding remarks. 



II. FLUCTUATION CONDUCTIVITY IN 2D 
SUPERCONDUCTORS 

The definition of fluctuation conductivity within the 
BKT theory relies on the Bardeen-Stephen formula 29 , 
which gives the excess conductivity Act = a — <jn with 
respect to the normal-state conductivity <jn as a function 
of the density of free vortices Hf above Tbkt as: 



Act 



(1) 



where the vortex mobility is [iy = 27r£ 2 c 2 p„/$ 2 l , £o is 
the spacing for the vortex lattice, that we will assume 
equal to the zero-temperature coherence length, p n is the 
normal-state resistivity and $o the flux quantum. The 
vortex density is also conventionally defined in terms of 
the correlation length £ as 2ixnp = l/£ 2 , so that near 
the transition where Act 3> on the ratio between the 
resistance R and its normal-state value Rn is: 



R 

R~ 



N 



(2) 



which is the formula usually quoted in the literature for 
the paraconductivity due to vortices, once that £(T) is 
calculated within BKT theory. Remarkably, in 2D the 
same formula is valid for the Aslamazov-Larkin contri- 
bution of GL fluctuations 16 ' 30 . However, in this case 
the temperature dependence of £(T) follows £q L ~ 
l/log(T/T c ), or equivalently £% L ~ T C /(T - T c ), where 
T c is the GL or mean-field critical temperature. In 2D 
films the true transition occurs at a Tbkt lower than 
T c , the distance between the two being a function of the 
microscopic parameters of the system. Thus, the fluctu- 
ation conductivity crosses over from a GL regime, where 
it shows a tendency to a power-law divergence at T c , to a 
KT regime, where the correlation length diverges asymp- 
totically when T — > Tbkt as £ ~ exp(b/\/i), where we 
introduced the reduced temperatures: 



t 



T — Tbkt 
Tbkt 



tc 



T c — Tbkt 
Tbkt 



(3) 



As it was observed already by Halperin and Nelson 17 , 
the BKT fluctuations are expected to be present only 
in the range of temperatures t <C t c . Thus, using the 
general expression Act oc £ 2 , they proposed a well-known 



interpolation formula between the GL and BKT regimes 
given by: 



Act = a hn cr N smh 2 {yjb hn t c /t) 



(4) 



where b^n is a dimensionless constant of order one. Thus, 
for bhntc/t S> 1 one recognizes the exponential divergence 
characteristic of the BKT theory, while for bhntc/t <C 1 
one recovers a power-law increase typical of GL fluc- 
tuations. For the prefactor a; m Halperin and Nelson 
assumed a^„ = 0.37/6; m , following the Beasley-Mooij- 
Orlando 31 approximate relation between on and the ra- 
tio Tbkt/T c in dirty BCS superconductors. We note 
in passing that in the dirty BCS limit one can indeed 
consider only the AL contribution to the GL paracon- 
ductivity, while in clean samples also the Maki-Thomson 
contribution can be sizable, making the analysis of GL 
fluctuations more involved 19 . 

In the HN formula some ambiguity is still present in 
the choice of the parameters, i.e. the prefactor and the 
exponential coefficient. In this work we want to fix this 
ambiguity by determining exactly their values from a RG 
analysis of the BKT correlation length. In this way, we 
provide a clear procedure to analyze experimental data 
by respecting the internal consistency between the fitting 
parameters for the paraconductivity. Moreover, we will 
discuss how finite-size effects and inhomogeneity affect 
the fluctuations resistivity above and below Tbkt- All 
these issues turn out to be crucial to correctly interpret 
the experimental data, as we will discuss in Sec. IV and 
V. 



III. BEHAVIOR OF THE CORRELATION 
LENGTH WITHIN THE BKT THEORY 

As it was shown by Kosterlitz in its original work 12,13 , 
the critical properties of the BKT transition can be cap- 
tured by the analysis of the RG equations for the two 
main quantities involved in the SC transition: the super- 
fluid stiffness J and the vortex fugacity g — 2ne^ l3 ^ T ' > 
(f3 = 1/ksT), where /x is the vortex-core energy. The 
stiffness J is the energy scale associated to the 2D su- 
perfluid density n 2D (measured experimentally via the 
inverse penetration depth A), including already the tem- 
perature depletion due to quasiparticle excitations: 



J(T) 



H 2 n 2D 



(T) 



Am* k 



B 



167re 2 A 2 (T)' 



(5) 



where d is the film thickness and m* is the effective mass 
of the carriers. The vortex-core energy fi is in general a 
multiple of the stiffness itself 21,24 ' 32 : 



H(T) = jlJ(T) 



(6) 



where ft is a dimensionless constant. As we discussed 
recently 21 ' 24 ' 25 , a general approach to the BKT transition 
can require to assume that p, deviates with respect to the 
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conventional value j2xY = 7r 2 /2 that it acquires in the 
XY model 12 ' 13,32 , so that we shall use 



/i = a/ixy = a- 



(7) 



The RG equations in the variables g and K — nJ s /T 
can be written (in analogy with the notation used for the 
sine-Gordon model 14 ) as: 



(8) 
(9) 



where I? = lna/£o, and £o, a are the original and rescaled 
RG lattice spacing, respectively. The physical value of 
the superfiuid density J s is determined by the limiting 
value of K under RG flow, i.e. J s = TK(£ -> oo)/n. 
In the low-temperature regime K > 2 so that the vortex 
fugacity scales to zero under RG flow and J s is finite, with 
a small renormalization with respect to the initial value. 
Instead at high temperature K < 2 the vortex fugacity 
becomes relevant, it diverges under RG flow and as a 
consequence J s scales to zero. The BKT temperature is 
the one where the above system of equations reaches the 
fixed point K = 2, g = 0, so that at Tbkt 



itJs{Tbkt) 



T 



BKT 



(10) 



i.e. one recovers the universal jump of the superfiuid 
density 13-15 . 

The usual definition 12 ' 13 of the correlation length £ 
above Tbkt relies on the determination in the RG equa- 
tions (8)- (9) of a characteristic scale £ at which the vor- 
tex fugacity is "sufficiently large". The exact definition of 
this scale is somehow arbitrary, but it does not change 
qualitatively the results. In practice, we shall use as a 
working definition the scale £ s where the RG parame- 
ter K{t s ) related to the superfiuid density vanishes. The 
vortex density and correlation length are then defined as: 



2imp 



9® 



1 



27ra 2 (^) £ 



72> «W 



doe* 



(11) 



where £ = rnin(£ s ,£ ma x), where £ m ax takes into account 
finite-size effects, as we shall discuss below. Near T B kt 
one expects to recover the well-known 12 exponential be- 
havior of £(T), that we will parameterize as: 



1 - L p b/Vt 

£o A 6 ' 



t -> 



(12) 



where b and A are constants of order one. 

It should be emphasized that if the bare superfiuid 
stiffness were a constant independent on the temperature 
Eq. (12) would be valid until t ~ 0(1), i.e. essentially at 
all relevant temperatures above Tbkt- However, what 
limits in a crucial way the applicability of the above ap- 
proximation is the temperature dependence of J(T) due 



to quasiparticle excitations, which lead to the vanishing 
of J(T) at the GL temperature T c . In a s-wave BCS 
superconductor a good approximation for J(T) is: 



AT) 
Jo 



m 



A(T) = A tanh 




For any practical purpose, to determine t c what matters 
is only the behavior of J{T) near T c , that according to 
the above equation is linear: 



J{T) « Jo T I I 



T 



T « T r . 



(13) 



By neglecting the normalization of J s with respect to J 
due to vortices already below Tbkt, one can approxi- 
mately estimate the BKT temperature by the condition 
(10) J {Tbkt) = 2T B kt/^, so that one gets: 



tc 



ST C 

7T 3 J ' 



(14) 



Thus, one sees that the larger is Jq/T c the smaller 
is the interval t c . In a thin film of conventional 
superconductors 18-20 the typical mean-field tempera- 
tures are of order of few K, while J can be as large 
as the Fermi energy if n 2D (0) coincides with the elec- 
tron density, as it is the case in clean superconductors 29 . 
Indeed, Eq. (5) gives: 



J = 22.1 n 2Z3 [10 13 cm -2 ] K = 0.62 



d[A] 



A 2 [(^m)] 



K (15) 



With d of order of few nanometers and A(0) <~ 100 A, as 
it is the case for conventional clean superconductors 29 , 
one has Jo of order of 10 3 K, so that t c would be of order 
of 10~ 3 , and then the BKT transition would be essen- 
tially indistinguishable from the mean-field T c . However, 
in dirty films of superconductors Jo can be substantially 
reduced with respect to the clean case 18-20 , so that t c as- 
sumes usually values around 0.05 < t c < 0.5, and Tbkt 
is sufficiently far from T c to be observable. 

From the point of view of the temperature dependence 
of the correlation length one sees that as one moves away 
from Tbkt toward T c the decrease of n(T) near T c makes 
the increase of the vortex fugacity under RG flow very 
fast, so that the RG superfiuid density vanishes for a{£) ~ 
£o and £ scales as the unrenormalized vortex fugacity 24 , 
i.e. 



-0fi(T)/2 



t -> t c 



(16) 



Since fi(T) -> as T -> T c , it follows that £ = Co at 
T = T c . Notice that strictly speaking one could as- 
sume that also the vortex-lattice spacing £o increases 
as T — > T c as the mean-field correlation length. How- 
ever, this just signals that as T c is approached the phase- 
modulus separation that justifies the BKT approach to 
the fluctuations fails, because the two degrees of freedom 
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cannot be no separated anymore in a controlled way. One 
expects that £(T) interpolates in a continuous way be- 
tween the BKT regime (12) and the GL regime, as we 
shall see below. Anyway, what is crucial to realize is 
that the approximate form (12) is limited to a regime 
t -C t c by the existence itself of a fluctuating GL regime 
at higher temperatures. 

As it is well know 12 , in an infinite system the BKT cor- 
relation length remains infinite anywhere below Tbkt- 
This signals the fact that in 2D no symmetry breaking 
can occur, so that the SC correlations decay to zero at 
large distance with a non-universal power-law decay, in- 
stead of the exponential decay to a finite value that one 
has in higher dimensions. From the point of view of the 
definition of £ used above in terms of the scale £ s where 
J s (£) vanishes, since J s is always finite below Tbkt then 
£ = oo in the SC phase. However, in any real system the 
RG scaling must be stopped at a certain scale £ max 



£max — kog * ) 
SO 



(17) 



where L is the maximum physical scale accessible in the 
system (like the system size or the size of the homoge- 
neous domains, see Sec. IV). As a consequence, the di- 
vergence itself of the correlation length is cut off below 
the temperature where £ s = £ max , so that the correla- 
tion length starts deviating from Eq. (12), and finite-size 
effects dominate. Sufficiently below Tbkt K(£) is sub- 
stantially unrenormalized with respect to the initial value 
ivf(O), 14 so that we can estimate finite-size effects by inte- 
grating the RG equation for g using K(£) ~ K(0). Thus, 
g u (£) — g u (0)e^ K ' e from Eq. (9) and the correlation 
length behaves as: 



Co 



= e MT)/2 



K(T)/2 



T < Tbkt- 



(18) 



We notice that, following Eq. (13), near T c , K(T) can be 
approximated as K(T) = 2(1 — t/t c ). As a consequence, 
from Eq. (18) we see that the size-limited £ grows below 
Tbkt as: 



1 

£o 



(19) 



Thus, if one considers two cases with the same L and 
Tbkt but different t c the finite-size effects are less pro- 
nounced in the case with smaller t c , because £ is still very 
large below Tbkt- 

As an example of the correlation-length temperature 
dependence we show in Fig. 2 the numerically-integrated 
£ at various values of the scale L, along with the fit of 
Eq. (12). Here we consider the case Jq = T c = 0.25 K, 
so that a relatively large interval t c ~ 0.36 is obtained, 
to better shown the BKT to GL crossover. From the 
inset of Fig. 2 we observe that to correctly capture the 
exponential scaling one must use extremely large sizes, 
far beyond the experimentally accessible regime. Indeed, 



since for conventional superconductors usually £o — 100 
A and L ~ 1 cm, then L — 10 6 £o and £ max — 12, while 
the exponential fit in Fig. 2 is better defined with £ max — 
100. 
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FIG. 2: (Color online) Temperature dependence of the cor- 
relation length within KT theory, obtained by numerical in- 
tegration of the RG equations (8)- (9) with the definition Eq. 
(11), where £ is determined either by the vanishing of the su- 
perfluid stiffness or by the finite system size f mal , see Eq. 
(17). Here T c = J = 0.25 K, and /} = jixv- The solid line 
represents the fit done using Eq. (12) with A=5 and b = 1.25. 
Since Tbkt = 0.18 the 6 value obtained by the fit is consistent 
with Eq. (21). Inset: same data plotted as a function of l/\/t, 
where t is the reduced temperature defined in Eq. (3). Notice 
that the curves for lower l m ax deviate from the infinite-length 
scale limit (12) at higher temperatures. 



As one can see in the main panel of Fig. 2, the expo- 
nential fit (12) deviates slightly from the RG correlation 
length as one moves to temperatures higher than Tbkt- 
This qualitative difference is also reflected by the tem- 
perature dependence of the fluctuations resistivity, that 
can be calculated according to Eq. (2), by using either 
the RG or the approximate form (12) of the correlation 
length. The result for the same set of parameters used in 
Fig. 2 is shown in Fig. 3: as shown in the inset, on a linear 
scale R has always an upward curvature, with no appre- 
ciable differences for the different L values. Such a differ- 
ence becomes instead evident on a logarithmic scale: as 
T — > Tbkt-, RJ Rn deviates from the infinite-size limit, 
exactly in the same fashion observed in experiments in 
thin films (see for example Fig. 9 of Ref. 19). This scaling 
can be probed by small magnetic fields: indeed, a small 
field has the effect to cut-off the RG at a scale Lb oc 1/B. 
The curves shown in Fig. 11 of Ref. 19 for the fluctuation 
resistivity at several (small) values of the magnetic field 
reproduce exactly the behavior observed in Fig. 3. 

Due to the need of exceedingly large system sizes for 
an accurate determination of the parameters of the fit 
(12), an a-priori estimate of them as a function of physi- 
cal parameters would be particularly useful. To address 
this issue we computed £ for several possible values of t c 
and //, to analyze the variations of b. While different t c 
values naturally occur in different systems, the choice of 
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FIG. 3: (Color online) Fluctuation resistivity in KT theory 
as a function of temperature for several values of the system 
size. The same set of parameters of Fig. 2 has been used. 
Inset: same curves but on a linear scale, where no appreciable 
difference can be observed for the different system sizes. 



fl is somehow still an open problem. As we mentioned 
above, within the XY model (that is one of the possible 
models for describing phase fluctuations in a supercon- 
ductor) jlxY = 7r 2 /2: this value naturally arises from 
the mapping of the discrete XY model on the continuum 
Coulomb-gas model 32 , and it takes into account the fluc- 
tuations at a scale of the order of the lattice spacing. In 
a BCS superconductor one could instead fix the value of 
the vortex-core energy by computing exactly the energy 
per unit- length of a vortex line 29 : 



( 


2 


[l A 






\ A " 


) 


log — - 

4o 


l-e 


= IT J 


log — + e 

4o 











so that according to our definition p, = 7re. A precise esti- 
mate of e ~ 0.497 for the vortex core in three-dimensional 
geometry is given in Ref. 33,34, so that within BCS one 
could eventually expect smaller values of ji, 



into account Eq. (12) and (16), we also see that the for- 
mula (12) can only be used when b/y/t 3> 1, which means 
t <C 4a 2 t c . Indeed, out of this regime £ decreases accord- 
ing to Eq. (16), and afterward (T > T c ) one enters the 
GL fluctuations regime. We notice also that this result 
agrees with the HN result (4), once one uses bhn = 2a. 




FIG. 4: (Color online) Estimate of b extracted by fitting the 
RG £(T) near Tbkt with Eq. (12) for several values of t c , a. 
The straight line is y = 2.1a;. 

The analysis we have made by coupling the standard 
BKT formula (12) with the RG analysis of the BKT tran- 
sition thus allows us to get strong constraints on what 
the "fit parameters" that are to be used in the BKT for- 
mula can be. This will help in the following to take into 
account additional effects such as the ones of inhomo- 
geneities, but specially to know if the fit to the BKT 
functional form corresponds to a physical fit, with rea- 
sonable parameters, or if one just "forces the fit" (see 
discussion in Sec. V). 



IV. ROLE OF INHOMOGENEITIES 



Hbcs — 2 



fJ-XY 



(20) 



Finally, we notice that recent analysis in the context 
of cuprate superconductor 21,24 ' 26,27 has explored instead 
the possibility that fi is larger than [Ixy ■ On the light of 
the previous observation, we considered the behavior of 
the correlation length for a range of values 0.5 < Ji < 1.2, 
and we extracted the corresponding b parameter by fit- 
ting data near Tbkt with Eq. (12). The results are 
shown in Fig. 4. As one can see, within a certain de- 
gree of uncertainty, b is found to scale approximately as: 



b ~ 2aVtc, 



(21) 



where a is the scale of the vortex-core energy defined 
in Eq. (7). Eq. (21) is the first important result of this 
paper. Indeed, it establishes a precise relation between 
the parameter b that appears in the typical exponential 
expression for the BKT correlation length and the dis- 
tance t c between the GL and BKT temperature. Taking 



Once we established a clear framework for taking 
into account the exponential behavior of the correlation 
length and the finite-size effects we can improve the HN 
original interpolation formula for paraconductivity, in or- 
der to obtain a self-consistent treatment of SC fluctua- 
tions all the way below and above Tbkt- We propose 
the following interpolating formula: 



R 



1 



1 



1 + (A<7/<7„) 

where £ is given by: 



i + (£/£o) 2 



(22) 



£ 2 6 
— = — sinn —p , 

£o A y/i 



■kJ(T)/2T 



T > Tbkt 



T<T B kt (23) 



(24) 



with the parameter A, b obtained by the calculated be- 
havior of the numerical RG correlation length near the 
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transition, so that b is given by Eq. (21), and A is a 
number of order 1. As it is shown in Fig. 5 this is in- 
deed a very good approximation for the numerical RG 
solution near and below Tbkt- Moreover, even though 
it does not capture the regime (16), where the estimate 
(24) is larger than the RG correlation length, it correctly 
reproduces the GL fluctuation regime at T » T Cl where 
Act ~ T C J (T — T c ). We stress once more that within this 
approach the temperature dependence of the correlation 
length in the whole fluctuating regime is uniquely de- 
termined by the critical temperatures Tbkt,T c and the 
value of the vortex-core energy, that can be fixed by the 
comparison with the experimental resistivity data. 




0.25 



FIG. 5: (Color online) Comparison between the numerical RG 
correlation length defined by Eq. (11) and the approximate 
formula (23)-(24). Here we used the same Tbkt ,T c , fj, values 
of Fig. 2 in the case imax = 5, so that b — 1.25 and A — 5. 

The main physical message of the set of Eqs. (21), (23)- 
(24) is that from the temperature dependence of the fluc- 
tuation resistivity between Tbkt and T » T c one can 
deduce interesting informations also on the microscopic 
parameters of the superconducting system, that deter- 
mine the distance between T c and Tbkt- To illustrate 
the application of this approach we consider the analysis 
of the resistivity data in superconducting heterostruc- 
tures reported in Ref. 4. Here the resistive transition 
occurs around 0.19 K, but with a relatively large tail 
with respect for example to what observed in 2D films 
of ordinary superconductors 18-20 . If we neglect the tail, 
by means of Eqs. (23)- (24) we can obtain the curve la- 
beled as 'Horn' in Fig. 6. As one can see, we reproduce 
the overall shape of the resistivity, except from the tail. 
Since the fit gives b = 0.19, according to Eq. (21) we de- 
duce that T c and Tbkt almost coincide, with T c = 0.19 
K and Tbkt = 0.188 K. 

Let us discuss now the origin of the remaining tail of 
the resistivity near the transition. First of all, we notice 
that it cannot be due to finite-size effects, that can be 
treated exactly within our approach. Indeed, even using 
a relatively small L = 2 /mi (as suggested by the critical 
current, see below), with £o = 70 A 4 , one gets £ max = 
5.6, but due to the small t c value the £ from Eq. (23) is 
still very large below Tbkt, so the finite-size effects by 



themselves are not responsible for the observed tail. On 
the other hand, from Fig. 3 we notice that a tail with 
upward curvature is typical of the fluctuation resistivity 
near the Tbkt transition. However since t c in this case 
is extremely small, the tail cannot be even resolved in the 
scale of Fig. 6. On the contrary, if the transition itself 
is broadened one could expect to enhance the BKT tail 
and to reproduce the experimental data. 

On the light of this observation, we suggest that the 
resistive tail can be attributed to an inhomogeneous spa- 
tial distribution of the local superfluid, to be ascribed to 
an intrinsic inhomogeneous density distribution in these 
systems. In analogy with the the analysis of inhomo- 
geneity on the superfluid density performed in Ref. 21, 
we shall assume for simplicity a gaussian profile: 



P(J) 



1 



2nS 



exp 



-{J -.If 



2S 2 



(25) 



where J is the Jo value determined by the above fit for the 
homogeneous case. To each J corresponds a mean-field 
temperature T rn f — T C (J/J), and accordingly a given 
Tbkt value. While far from Tbkt such an inhomogene- 
ity is harmless, near Tbkt it will give important effects, 
once we average over different patches having different 
transition temperatures. To compute the effect of the 
sample inhomogeneity on the resistivity we go back to 
Eq. (2), that defines the resistivity due to vortices. We 
notice that the main quantity that enters the theory is the 
vortex density: indeed, the correlation length defined by 
Eq. (2) is just a different way to express the vortex den- 
sity. Thus, in the presence of the inhomogeneity (25), we 
must average the vortex-density values np{J) obtained 
in each patch with a given local J value: in this way, even 
below Tbkt there will be patches of the system where 
up is finite, because J < J and consequently the local 
Tbkt is smaller than the average one. As a consequence, 
the BKT tail gets enhanced, in excellent agreement with 
the experiments (with S = 0.02J), see Fig. 6. 

A second effect of the inhomogeneity that emerges in 
the experiments of Ref. 4 is the lack of the universal jump 
of the superfluid density near J s . In Ref. 4 J s has been 
indirectly measured via the I — V characteristics in the 
non-homic regime (see discussion above Eq. (28)), given 
within the KT theory by 13 - 17 ' 35 : 



y = T l+-K.J S (T)/T _ ja{T)_ 



(26) 



Thanks to the universal relation (10), right below the 
transition ttJ s (Tbkt)/Tbkt = 2, so that the coefficient 
a = 3. On the other hand, in the inifite-size and homo- 
geneous case J S (T) jumps discontinuously to zero right 
above Tbkt, so that the coefficient a(T) is expected to 
jump discontinuously to 1. Even though finite-size effects 
remove the jump and lead to a rapid downturn (see inset 
of Fig. 6), the measured temperature dependence of a(T) 
is still much smoother. This effect can be explained by 
the presence of inhomogeneity, that leads to a smooth 
downturn of J s near Tbkt 21 - 




0.16 0.18 0.2 0.22 0.24 0.26 0.28 0.3 
T(K) 



FIG. 6: (Color online) Comparison between the resistivity of 
the heterostructure measured in Ref. [4] and the resistivity 
obtained with the interpolating formula (22). The curve la- 
beled 'Horn' refers to the case of a single local J value, while 
the curve 'Inhom' refers to the resistivity obtained by sample 
average over the local distribution (25) of superfluid density. 
Inset: coefficient a(T) = I + tt.J s (T)/T of the J - V charac- 
teristic in the two cases. 



T(K) 




FIG. 7: (Color online) Comparison between the experimen- 
tal data in the unbiased sample (V=0) of Ref. [5] and the 
resistivity obtained with the interpolating formula (22), with 
('Inhom') and without ('Horn') inhomogeneity. The parame- 
ter values are given in the text. The two arrows mark Tbkt 
and T c obtained with the fit in the homogeneous case. 



As a further check of the correctness of our analy- 
sis of the data, we discuss the implication of the in- 
terpolated GL+BKT fit on the estimated value of the 
superfiuid-density value at T = 0. From the value of 
t c « 8T c /(tt 3 J ) = 0.0106 obtained by the fit, it fol- 
lows (see Eq. (14)) that J = 7 K. According to Eq. (5) 
above, Jo is in turn controlled by the zero-temperature 
value of the density of superfluid electrons. While in a 
clean superconductor n 2D = n 2D at T = 0, where n 2D 
is the sheet carrier density of the sample, in the dirty 
case only a fraction of the order of £d/£,o of the electron 
density condenses into the superfluid fraction at T = 0, 
where l& = vft is the mean-free path, with vf Fermi ve- 
locity and r scattering time. Since the BCS correlation 
length at T — is £o = Hv p j 'tc A(0) , with A supercon- 
ducting gap, and r can be determined by the sheet resis- 
tance Rq = n 2D e 2 T/m* (here □ denotes the sheet area), 
we obtain an estimate of the disorder-reduced superfluid 
density as: 

h 2 n 2 s D _ h 2 n 2D ttA(0)1 _ 
B ° 4m* 4m* Hvf 

H n 2D e 2 T ttA(0) _ Re ttA(0) 

m* 4~ ~ R^ 4 ^ ' 

Using the universal value R c = h/e 2 — 41140/D, the 
measured normal-state resistance R ~ 390O/D, and the 
BCS estimate A(0) = 1.76T C with T c = 0.19 K, we obtain 
Jo ~ 3 K, which is in excellent agreement with the value 
J = 7 K obtained by the fit, taking into account the 
approximate estimate of the condensate fraction used in 
Eq. (27) above. 

The crucial role of inhomogeneities in controlling the 
tails of the resistive transition is made even more evi- 
dent by extending the above analysis to a second example 



of superconducting heterostructures measured in Ref. 5, 
where the normal-state resistance is higher, so that inho- 
mogeneous effects could be expected to be stronger. We 
shall consider the case of zero gate voltage in Ref. 5, i.e. 
the as grown sample without the field-induced doping. 
The data are shown in Fig. 7 along with our proposed fit 
(T c = 0.269 K, T BKT = 0.25 K, b = 0.63 and J = 1.05 
K). Notice that a larger b is consistent with a larger t c , 
that in turn can be attributed to a larger disorder in this 
second case. Indeed, in this sample Rq = 2400 fi/CJ, 
and since n 2D — 4.5 x 10 13 cm -2 using Eq. (27) we ob- 
tain a consistently smaller estimate of Jo ~ 0.6 K, which 
is exactly the result of the fit (a factor 10 smaller than 
the previous case). Moreover, the inhomogeneity is also 
larger, with 6 = 0.1 4 J, and this explains the very large 
tail of the transition. 




FIG. 8: (Color online) Schematic sketch of the mesoscopic 
structure of a 2D film that can account for the inhomogeneity 
effects discussed in this work. Even though the system does 
not have a true granular structure the homogeneous regions 
will have a typical size Lhom smaller than the physical size 
iphys of the sample. This explains why L = Lh om < L p h ys 
in Eq. (29), giving rise to a critical current for linear-to- non- 
linear characteristic larger than expected for the homogeneous 
case. 
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As we emphasized above, finite-size effects alone can- 
not account for the broadening of the transition. How- 
ever, one could convert in a hand- waving way the in- 
homogeneities effects that we discuss into some form of 
typical size for "homogeneous" grains in the system. Let 
us emphasize that this is just a way to visualize the role of 
the inhomogeneities and that we do not assume any sharp 
granular structure here, see Fig. 8. The inhomogeneity 
distribution (25), and in particular the broadening of J 
values, (5, could be converted into a characteristic size 
ihom much smaller than the true system size. An indica- 
tion in this sense is provided by the analysis of the I — V 
characteristics reported in Ref. 4 . Indeed, within BKT 
physics the I —V characteristic is ohmic at low current 
and non-homic when the external current is sufficiently 
high to dissociate vortex-antivortex pairs. Below Tbkt, 
the critical current I* above which current-induced free 
vortices lead to the anomalous power-law dependence de- 
scribed by Eq. (26) above, is in turn related to a char- 
acteristic length scale L where finite-size effects become 
relevant 13,17,35 , so that: 

r = 27r j s ^^ (28) 

where L p h ys is the physical dimension of the sample. By 
using the approximate relation between J s and Tbkt, 
one can estimate for the sample of Ref. 4 that 

/* ~4K B T BKT ^^^ = ^^0.5xlO- 8 A (29) 
$o L L 

Since experimentally the critical current is of order of 
0.5 x 10~ 6 A and L p hy S — 0.2 mm, we deduce that 
L ~ 2/xm. This can be in turn identified with the size 
Li lom of the homogeneous domains, so that £ max = 5.6, 
as we used in our calculations. Notice that Eq. (28) ac- 
counts also for the temperature dependence of I* (T) ob- 
served in the experiments: indeed, as we have seen J S (T) 
increases rapidly below Tbkt (see the behavior of the 
a(T) exponent in Fig. 5), leading to a very fast increase 
of the critical current below the Tbkt transition. At 
temperatures well below Tbkt one should further ac- 
count in the analysis of the I — V characteristic for the 
nucleation of virtual vortex-antivortex pairs, formed by 
a single vortex and its virtual image due to the edge cur- 
rents in a finite-size system, as pointed out recently in 
Ref. [36]. 

V. DISCUSSION ON THE ANALYSIS OF THE 
EXPERIMENTAL DATA 

In the previous Section we performed a systematic 
analysis of the resistivity data from Refs. 4,5. Our analy- 
sis is based on the regular BKT transition and its connec- 
tion with Gaussian fluctuations leads to the conclusion 
that the BKT transition occurs very near (within few 
percent) of the mean- field GL transition temperature T c . 



In this case the BKT fluctuation regime is restricted to a 
small range t c of reduced temperatures near Tbkt, where 
the correlation length (24) has the exponential behavior 
(12), with an exponential divergence near Tbkt that is 
controlled by the same parameter b (21) that measures 
the distance from the GL transition T c . At temperatures 
near and below Tbkt inhomogeneities are responsible for 
the tail in the resistive transition. 





Data 










Fit Eq. (30) - 






















/■' 1 


















' 0.8 
0.6 
0.4 


/ b"=i — ; 

; b=o.i -- zz ^_- 








0.2 





















3 | 1 2 3 
~t\ex t 


4 



0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
T(K) 

FIG. 9: (Color online) Example of application of the the BKT 
fit proposed in Ref. [28] for the same data shown in Fig. 7. 
The theoretical curve is Eq. (30) with Tbkt = 0.245 K and 
b — 0.106. Inset: temperature dependence of the function 
f(t) = exp(— 2b/ yt) for two b values. For 6=1 one can 
clearly distinguish on this scale a first regime with upward 
curvature, up to the temperature Tfi ex given by Eq. (31), 
followed by a regime with downward curvature. For b = 0.1 
only the second regime is visible on this scale. 

A quite different way to interpret the data was followed 
instead in Refs. 4 ' 5 ' 28 , based essentially on the idea that 
T c is far larger than Tbkt, so that the whole fluctua- 
tion regime should be described via the standard BKT 
approach, see Fig. 1. For the sake of clarity in the fol- 
lowing we will refer to our approach as 'BKT-GL'(Fig. 
la), and to the one proposed in Refs. 4 ' 5 ' 28 as 'BKT-only' 
(Fig. lb). In the latter case, the expression (12) for the 
correlation length should be valid in all the temperature 
regime where R/Rn deviates from 1, and b is a number 
of order 1, to assure that T c is enough far from Tbkt- 
As a consequence, the whole resistivity (2) is given at all 
temperature above Tbkt by: 



-5- = A 2 e- 2b V^KT/(T-T BKT ) = A 2 (30) 
Rn 

where f(t) = exp(— 2b/y/t) and t is the reduced temper- 
ature (3). An example of application of Eq. (30) to the 
same resistivity data analyzed in Fig. 7 above is presented 
in Fig. 9. Here we used Eq. (30) with Tbkt = 0.245 K 
and considered for the moment b as a completely free pa- 
rameter. In this reasonable fit of the data can be 
obtained by using b = 0.106. Such a small value of b can 
be understood by looking at the behavior of the function 
f(t) shown in the inset of Fig. 9. As one can notice, if one 
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extend f(t) at arbitrary t values above Tbkt it actually 
saturates to a constant value. The crossover from the 
low T regime where f(t) displays an upward curvature 
to the one where it displays an downward curvature can 
be estimated by the flex of this curve, that is found at: 

T fiex =T BKT [l + 0A4b 2 ]. (31) 

As a consequence, if b is very small the function f(t) sat- 
urates very rapidly and one can fit the experimental data 
with the expression (30) by assuming a very small b value. 
However, this procedure is totally inconsistent with the 
fact that within BKT theory b is not a free parameter, 
as we discussed in the previous Sections, but it depends, 
through Eq. (21), on the distance t c between the BKT 
and GL temperatures. In particular, the fit presented in 
Fig. 9 gives b = 0.1, that would imply a t c — 0.0025, in 
clear contradiction with the a-priori assumption that T c 
is so much larger than Tbkt than the whole fluctuation 
regime has BKT character. Moreover, a smaller b value 
would also imply that J is very large (J ~ 26 K), that 
is again inconsistent with the estimate (27) based on the 
normal-state resistivity value (i.e. J <~ 7 K). We note in 
particular that this procedure leads to a value for the pa- 
rameter b that is about one order of magnitude smaller 
than what we obtained in our case (see Fig. 7), where 
b = 0.6. Indeed, in the BKT-GL approach most of the 
fluctuation regime is accounted by GL fluctuations, so 
that Eq. (30) is actually valid only in the limited range 
of temperatures T -C Tfi ex where f(t) displays an up- 
ward curvature. 

An incorrect use of the BKT formulas can lead to quite 
unphysical values for the b parameter, as it is the case in 
our opinion for Ref. [28]. Moreover, independently of the 
question of the distance between Tbkt and T c , a second 
drawback of this BKT-only analysis is that it attributes 
the tail of the transition, that is the real signature of 
BKT physics, to unrealistic finite-size effect. Indeed, in 
Ref. [28] it has been proposed that the deviation from the 
best fit obtained with Eq. (30) and the experimental data 
occurs when £ is cut-off by finite-size effects. As we have 
seen above, finite-size effects become relevant near Tbkt 
when £ ~ L, see Eq. (18). Thus, since R/R N = (£ /£) 2 
according to Eq. (2), in Ref. [28] it has been argued that 
the fit with Eq. (30) deviates from the experimental data 
when £(T) = L, with T > T B kt- By applying this idea, 
since from Fig. 9 it follows that the fit deviates from the 
data already at R/ Rn ~ 0.5, one would obtain L ~ 2£o ~ 
100 A. This is an unrealistic small number for the size of 
the homogeneous domains, because it would not allow at 
all the formation of coherent vortex structures (and then 
the observation of their signatures) on a distance as small 
as only twice the lattice spacing for phase fluctuations 
(i.e. £o)- In our approach instead the deviation of the 
experimental data from the homogeneous-case fit is due 
to inhomogeneity, while finite-size effects alone cannot 
account for it. 

The only way to possibly justify the BKT-only analysis 
of the data would be if for some reason the constraints on 



the b parameter that we establish in our paper for a con- 
ventional BKT transition would not be obeyed, namely 
if the transition is of a different nature than the con- 
ventional BKT transition. What could be such an al- 
ternative is unclear from a theoretical point of view. A 
proposal in this sense has been formulated in Ref. [4], 
where it has been argued that in these systems one can 
estimate the initial value of the vortex fugacity in the 
RG equations (8)- (9) to be quite large, invalidating the 
applicability itself of the perturbative RG 13,37 . In this 
case, following the analysis of Ref. [37], the BKT transi- 
tion has a qualitative different character, since it refers to 
the melting of the 2D lattice of vortices formed at high 
vortex density. Such a scenario can provide for exam- 
ple an alternative interpretation to the high values of the 
critical current I* (28) mentioned above (see Supplemen- 
tary Information of Ref. [4]). This is thus an interesting 
proposal, but it leads to a BKT-only analysis of the data 
that still suffers from two main problems: a) No theoret- 
ical work has been devoted so far to explore the signa- 
tures of such a melting transition, if it exists, on the SC 
fluctuations contribution to the resistivity. In particular 
it is not at all obvious that such a scenario would lead 
to a different b parameter, so that the small b value ob- 
tained by the application of Eq. (30) at all temperatures 
above Tbkt would be meaningful, b) Recent numeri- 
cal simulations on the vortex liquid do not confirm the 
existence of this crystal phase 38 at even at high vortex fu- 
gacity and only standard BKT vortex phase seems to be 
present, so that it is not even clear whether such a melt- 
ing transition would exist. Thus, more theoretical and 
experimental work is certainly needed to test this possi- 
bility and to decide which one of the interpretations, the 
BKT-GL approach or the BKT-only based on the idea 
of vortex lattice melting, is indeed the correct explana- 
tion of the existing data. For example, further measure- 
ments of the size-dependence of the critical current I* 
or the direct measurement (via two-coils mutual induc- 
tance experiments) of the temperature dependence of the 
superfluid density can provide us with additional infor- 
mation on the system, to be tested against the different 
theoretical proposals. 



VI. CONCLUSIONS 

In the present work we investigated the role of finite- 
size effects and inhomogeneity of the BKT transition in 
quasi-2D SC systems. By analytical and numerical analy- 
sis of the RG equations for the BKT transitions we deter- 
mined in a precise way how the properties of the system 
control the behavior of the correlation length above and 
below the transition. By interpolating the BKT behavior 
with the GL fluctuations we propose an unified scheme 
to determine both the mean-field transition temperature 
T c and the BKT one Tbkt from a fit of the resistivity 
data. In our approach, those that are usually treated as 
free parameters in the BKT fitting formula are implicitly 
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determined by the these temperature scales, and allows 
us to identify the effects that must be attributed to the 
finite system size and inhomogeneity. The direct com- 
parison with recent experiments in SC heterostructures 
suggests that spatial inhomogeneity, whose role we dis- 
cussed already in the context of ultra-thin films of high-T c 
superconductors 21 , is a common ingredient to these sys- 
tems as well. Our analysis allows us to estimate also the 
superfluid-density content of these unconventional SC in- 
terfaces, whose direct experimental determination is not 
yet available. Interestingly, in Ref. [5] it has been shown 
that by decreasing the carrier density by field-effect one 
can induce in such heterostructures the transition from 
a SC to a metallic and then insulating state. An in- 
teresting open issue is the role played by inhomogeneity 
in such a crossover: for example, one could expect that 
as Jo(T c ) decreases and 5 increases one reaches the per- 
colative threshold below which the system cannot sustain 
anymore superfluid currents, in analogy with the behav- 



ior of diluted XY models 39 . This issue can have profound 
consequences on the nature of the quantum critical point 
inferred in Ref. [5], and certainly deserves further inves- 
tigation. 
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